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ABSTRACT 



, ^ We report on a sensitive search for H2 1-0 S(l), 1-0 S(0) and 2-1 S(l) ro-vibrational emission at 2.12, 2.22 and 2.25 fim in a sample 

' of 15 Herbig Ae/Be stars employing CRIRES, the ESO-VUT near-infrared high-resolution spectrograph, at R~90,000. We report the 

T-H detection of the H2 1-0 S(l) line toward HD 100546 and HD 97048. In the other 13 targets, the hne is not detected. The H2 1-0 S(0) 

T-H and 2-1 S(l) lines are undetected in all sources. These observations are the first detection of near-IR H2 emission in HD 100546. The 

H2 1-0 S(l) lines observed in HD 100546 and HD 97048 are observed at a velocity consistent with the rest velocity of both stars, 

I I suggesting that they are produced in the circumstellar disk. In HD 97048, the emission is spatially resolved and it is observed to 

P^ extend at least up to 200 AU from the star. We report an increase of one order of magnitude in the H2 1-0 S(l) line flux with respect to 

^^ previous measurements taken in 2003 for this star, which suggests line variability. In HD 100546 the emission is tentatively spatially 

resolved and may extend at least up to 50 AU from the star. Modeling of the H2 1-0 S(l) line profiles and their spatial extent with flat 

keplerian disks shows that most of the emission is produced at a radius larger than 5 AU. Upper limits to the H2 1-0 S(0)/ 1-0 S(l) and 

H2 2-1 S(l)/l-0 S(l) line ratios in HD 97048 are consistent with H2 gas at T>2000 K and suggest that the emission observed may be 

Q produced by X-ray excitation. The upper limits for the line ratios for HD 100546 are inconclusive. Because the H2 emission is located 

;_( at large radii, for both sources a thermal emission scenario (i.e., gas heated by collisions with dust) is implausible. We argue that the 

"^ observation of H2 emission at large radii may be indicative of an extended disk atmosphere at radii > 5 AU. This may be explained 

jrt by a hydrostatic disk in which gas and dust are thermally decoupled or by a disk wind caused by photoevaporation. 

Key words, circumstellar matter, stars - pre-main-sequence, emission-line, protoplanetary disks 
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transitions are very weak. In consequence, in contrast to other 



^ Circumstellar disks around pre-main sequence stars are the birth g^^ ^^^^^^^ (j^ost notably CO), H2 emission is harder to detect. 

^ places of planets. Hence, the characterization of their physical ^ith the advent of space borne observatories in the UV and in- 

1^ properties is of paramount importance. At the time when giant f^^^^^^ ^^^ g^Q^jn^ ^ased high-resolution infrared spectrographs, 

O planets are in formation, protoplanetary disks are rich in gas. j^ ^as been possible to start the search and study of H2 emission 

^ However, there are relatively few observational constraints of the f^om disks around young stars from the U V to the mid-IR. 
'""I gas content in the disk, in particular for the inner disk (R < 10 

^ AU), the region where planets are expected to form. To study In the UV, H2 electronic transitions trace in emission hot H2 

• ^H the gas content in the disk, molecular and atomic line emission gas in the disk that is photoexcited ("pumped") by Lya pho- 

/\^ is used. Because the disk has a radial temperature gradient, dif- tons (e.g. Valenti et al. 1993; Ardila et al. 2002; Herczeg et al. 

J-H ferent transitions of diverse gas tracers probe different radii in 2006) or excited by electrons that are generated by X-rays (e.g., 

^ the disk (see reviews by Najita et al. 2007 and Carmona 2010). Ingleby et al. 2009, France et al. 2011). In absorption, UV tran- 

These emission lines are in general produced in the surface lay- sitions trace cold, warm, and hot H2 in the line of sight (e.g., 

ers of the disk where the dust is optically thin. Martin-Zaidi et al. 2005, 2008a), either from the disk itself (e.g., 

By far the main constituent of the gas in protoplanetary disks flared disk, close to edge-on disk, disk wind) or the star's circum- 

is molecular hydrogen (H2). However, given its physical nature stellar environment (e.g., envelopes). In the mid-IR, H2 pure- 
rotational transitions trace thermal emission of warm gas at few 

5eW offprint «^t.v to: A. Carmona, e-mail: hundred K (e.g. Thi et al. 2001; Bitner et al. 2007; Martin-Zaidi 

andres.carmonaOunige.ch et al. 2007, 2008b, 2010; Lahuis et al. 2007; Carmona et al. 

* Based on observations collected at the European Southern 2008b). In the specific case of the near-infrared, ro-vibrational 

Observatory, Paranal, Chile. (Program IDs 079.C-0860C, 080.C- transitions of H2 trace thermal emission of hot H2 at a thou- 

0738A, 081.C-0833A) sand K, or H2 gas excited by UV or X-rays (e.g. Weintraub et 

** Postdoctoral Fellow of the Fund for Scientific Research, Flanders. al. 2000; Bary et al. 2003, 2008; Itoh et al. 2003; Ramsay Howat 



A. Carmona et al.: H2 near-IR emission from Herbig Ae/Be stars. 
Table 1. Coordinates, spectral types, distances, radial velocities (RV), and SED group " of the program stars. 



Star 



HD 58647 * 
HD 87643 '' 
HD 95881 
HD 97048 
HD 100546 
HD 101412 
HD 135344B 
HD 141569 
HD 144432 
HD 150193 
510ph 
HD 169142 
RCrA 
HD 179218 
HD 190073 



Q'(J2000.0) (5(J2000.0) Sp. Type rf [pc] RV[kms-'] Group" 

r^rr ^^ n^ , ^ , .. in a r^ ^ TAn-r-.T A ^ r. nW 77 



07:25:56.10 
10:04:30.29 
11:01:57.62 
11:08:03.32 
11:33:25.44 
11:39:44.46 
15:15:48.43 
15:49:57.75 
16:06:57.96 
16:40:17.92 
17:31:24.95 
18:24:29.78 
19:01:53.65 
19:11:11.25 
20:03:02.51 



-14:10:43.5 
-58:39:52.1 
-71:30:48.3 
-77:39:17.5 
-70:11:41.2 
-60:10:27.7 
-37:09:16.0 
-03:55:16.4 
-27:43:09.8 
-23:53:45.2 
-23:57:45.5 
-29:46:49.4 
-36:57:07.6 
+15:47:15.6 
+05:44:16.7 



B9IVe " 

B2e'' 
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B9Vne '' 

B9.5Ve ' 

F4Ve ■'■ 

AOVe'' 
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AlVe'' 
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A5Ve I' 

Al-F7e^ 

AOIVe '' 

A2IVpe ' 



28(r 

118^ 
180'' 
103'' 
160 « 

140^ 
99 rf 

145^ 
150'' 
130'' 
145^ 
130^ 
240'' 
>290s 



+36±2« 
+ 18±3'' 
+ 16±2'' 

-3±2« 
+2.5±1.5' 

-2±2 ^ 

-6±2« 
-11±3^' 

-3±2« 

0±2''' 

-9±2« 

+3±2s 



II 

I 

II 

I 

I 

I/II 

II, TD 

II, TD 

II 

n 

II 

I 
II 

I 
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Notes. '''' SED classification by Meeus et al. (2001). Group I sources are with "flared" disks and Group II sources are "self-shadowed" disks. 
TD means transition disk. These stars lack or have a weak near-infrared excess. This is interpreted as evidence of a hole or gap in the inner disk. 
Because the star shows evidence of accretion the gap refers only to the lack of small dust particles in the inner disk. **' Manoj et al. (2002) argue 
that HD 58647 may not be a Herbig Ae/Be star but a classical B[e] star. Baines et al. (2006) suggest that HD 58647 is a binary star based on 
spectropolarimetry of the Ha line. *''' The status of HD 87643 as a Herbig Ae/Be star is controversial. Oudmaijer et al. (1998) argue that HD 87643 
is an evolved B[e] star (see also Kraus 2009). Millour et al. (2009) based on AMBER- VLTI observations reported a close companion at 34 mas. 
References for the spectral type, distance, and heliocentric radial velocity: *'" van den Ancker et al. (1998); *''' SIMBAD; '-" Acke et al. (2004) and 
references there in; *^* Acke et al. (2005) and references there in; *'"' this work, see discussion Sect. 2.2; ''' Miiller et al. (2011); '-'* Kharchenko et 
al. (2007). 



& Greaves 2007; Carmona et al. 2007, 2008a). Because gas at 
these temperatures is expected to be located at radii up to a few 
AU, H2 near-IR emission has the potential of tracing the gas in 
the terrestrial planet region of disks if the observed emission is 
thermal. 

H2 near-IR emission from disks has been mainly studied 
toward the low-mass T Tauri stars, where it has been detected 
preferentially in objects exhibiting signatures of a high accretion 
rate (i.e., large Wa equivalent widths and U -Y excess) and in 
a few weak-line T Tauri stars with bright X-ray emission (see 
the discussion sections of Carmona et al. 2007, 2008a, and Bary 
et al. 2008). However, for intermediate mass young-stars, i.e., 
Herbig Ae/Be stars, there is to date only one reported detection 
of near-infrared H2 emission at the velocity of the star, namely 
in the star HD 97048 (Bary et al. 2008). Because Herbig Ae/Be 
stars are bright in the near-IR, the detection of the faint H2 
lines is challenging, in particular at low spectral resolution. 
In this paper, we present the results of a considerable effort 
that we undertook to search for H2 1-0 S(l), 1-0 S(0) and 2-1 
S(l) emission toward Herbig Ae/Be stars employing CRJRES, 
the ESO-VLT near-infrared high-resolution spectrograph, at 
resolution R ~ 90,000. 

The paper is organized as follows. In Section 2 we describe 
our CRJRES observations and the data reduction techniques. In 
Section 3 we present the resulting spectra and constrain the ex- 
citation mechanism of the detected lines based on H2 line ratios. 
We then analyze the observed H2 lines in the context of a flat disk 
model. In Section 4 we discuss our results. First, we compare the 
H2 line profiles detected with the line profiles of other disk gas 
tracers. We then discuss our H2 observations and observations of 
dust in the inner disk. Thereafter, we analyze H2 near-IR emis- 
sion in the frame of diverse disk structure models. We propose 
explanations for the relatively large number of non-detections, 
and explore diverse scenarios to explain the detections of near- 



IR H2 lines in HD 100546 and HD 97048. Finally, in Section 5 
we present our conclusions. 



2. Observations and data reduction 

2.1. Observations 

We observed a sample of 15 nearby Herbig Ae/Be stars of di- 
verse spectral types and disk geometries based on the spectral 
energy distribution (SED) classification introduced by Meeus et 
al. (2001) and SED modehng by Dullemond et al. (2001) (i.e., 
group I: flared disks, group II: self-shadowed disks). We sum- 
marize the properties of the stars of our program in Table [T] 

The sources were observed with the CRyogenic high- 
resolution InfraRed Echelle Spectrograph (CRIRES) mounted at 
the ESO-VLT UTl (Antu) atop Cerro Paranal, Chile (Kiiufl et 
al. 2004). The observations are standard long-slit spectroscopy 
observations. We employed a slit of width 0.2", with the slit 
rotated along the parallactic angle resulting in a spectral res- 
olution R of 90,000 (or 3.2 km s"'). We used adaptive optics 
(MACAO - Multi-Applications Curvature Adaptive Optics) to 
optimize the signal-to-noise ratio and the spatial resolution. The 
typical spatial resolution achieved was 0.22" (i.e., FWHM in the 
continuum). To correct for the sky emission, we used the stan- 
dard nodding technique along the slit employing a nod throw 
of 10". To correct for telluric absorption and flux-calibrate the 
spectrum, we observed spectrophotometric standard stars imme- 
diately before or after the science observations with the same in- 
strument setup. Our observations were carried out in service and 
visitor mode in 2008 under the ESO programs 080.C-0738A and 
081.C-0833A. We complemented our data for HD 97048 with 
CRIRES data taken by one of us with a similar setup (ESO pro- 
gram 079.C-0860C, PI E. Pantin). In these observations the slit 
was oriented in the north-south direction and we used a stan- 
dard star observed in our program 080.C-0738A to perform the 
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telluric correction and flux calibration. In Table |ATT| (in the ap- 
pendix) we present a summary of the observations and the point- 
spread function (PSF) FWHM achieved in the target and the cal- 
ibrator. 



2.2. Data reduction and calibration 

The data were reduced using the CRIRES pipeline V1.7.cQup to 
extraction of the ID spectrum, then custom IDL routines were 
employed for the telluric correction, accurate wavelength cal- 
ibration, barycentric correction, cosmic ray cleaning and flux 
calibration. For each chip flat-fielded image pairs in the nod- 
ding sequence (AB) were subtracted and averaged producing a 
combined image frame, thereby performing the sky-background 
correction. Raw frames at each nodding position were corrected 
for random jittering, employing the information in the fits head- 
ers. The ensemble of combined frames were stacked in one sin- 
gle 2D image spectrum, from which a one-dimensional spectrum 
was extracted by summing the pixels in the spatial direction in- 
side the PSF after background subtraction. The CRIRES pipeline 
provides a first wavelength calibration by cross-correlation with 
the Th-Ar lamp frame taken during the same night of the obser- 
vations, and the final product of the CRIRES pipeline is a ID 
spectrum. 

To coiTect for telluric absorption, the ID extracted science 
spectrum was divided by the ID extracted spectrum of the stan- 
dard star. The standard star spectrum was first corrected for dif- 
ferences in air-mass and air-pressure with respect to the science 
target spectrum employing 



hTD,„„.,„,j -k exp^- 



_^ ^TARGET PtARCET \ 
XsTD PsTD I 



Here X is the average of the airmass and P is the average of 
the air pressure. The continuum /q has been defined as uniform 
with a value equal to the mean plus three standard deviations 
(/o = /obs + 30"/^^,). The optical depth t is derived from the mea- 
sured data using t - -\vi{Iohslh)- There is some freedom in 
the choice of the continuum level, we selected the mean plus 
3cr to set the continuum over the noise level, which resulted in 
I obs I h) < 1 along the spectrum. The result of this connection is to 
make the depth of the telluric lines in the standard star spectrum 
similar to the depth of the telluric lines in the science spectrum. 
Small offsets of a fraction of a pixel in the wavelength direction 
were applied to the standard star spectrum until the best telluric 
correction (i.e., signal-to-noise) in the coiTected science spec- 
trum was obtained. 

By comparing the sky absorption lines in the non-telluric 
corrected ID science spectrum produced by the CRIRES pipe- 
line and a HITRAN model of Paranal's atmosphere, a final wave- 
length calibration was established. The typical wavelength ac- 
curacy achieved is 0.5 km s"'. The spectra were coiTected for 
the radial velocity (RV) of the star and the motion of the Earth- 
Moon-Sun system at the moment of the observation using the he- 
liocentric velocity correction given by the IRAPrltask rvcorrect 
and heliocentric radial velocities from the literature (see Table 

For HD 100546 and HD 97048, two measurements of the 
radial velocity are available from the literature. Donati et al. 
(1997) estimate an RV of 17+5 km s"' for HD 100546. Acke 
et al. (2005) derived an RV of 18 km s' for HD 100546, and of 
21 km s ' for HD 97048. To obtain an independent estimate of 

' http://www.eso.org/sci/data-processing/software/pipelines/index.html 
^ http://iraf.noao.edu/ 



the radial velocity of the objects, we used two additional meth- 
ods. First, we used archival high-resolution FEROSrl(R~40000) 
spectra of HD 100546 and HD 97048. We determined the ra- 
dial velocity by fitting Gaussians to narrow photospheric ab- 
sorption lines not affected by blending (because v sin /) or emis- 
sion, and comparing their centers to the center of Gaussians 
fitted to the same spectral features in rotationally broadened 
BLUERED (Bertone et al. 2008) high-resolution (R~500,000) 
synthetic spectral models corresponding to the spectral types of 
HD 97048 and HD 100546. We obtained for HD 100546 an RV 
of 16+2 km s', and for HD 97048 an RV of 17+2.5 km s^^. 
Second, we derived the radial velocity from the center of the CO 
ro-vibrational lines observed at 4.7 ^m in CRIRES R~ 90000 
spectra taken by our team (van der Plas et al. 2009). For this set 
of data we found an RV of 14±2 km s ' for HD 100546 and an 
RV of 16±2 km s' for HD 97048. 

Our results are consistent with the values of Donati et al. 
(1997) and the values of Acke et al. (2005), assuming a typical 
error of 2 km s"^ for the latter. Donati et al. (1997) derived their 
RV estimate for HD 100546 based only on the Mg 11 doublet at 
4481 A. Because we used 12 additional photospheric lines in our 
determination of the RV for HD 100546, we consider our estima- 
tion to be more precise (indeed, we obtained an RV of 17.3 km 
s ' for the Mg 11 lines). The principal limitation on the derivation 
of the radial velocities using optical spectra is that we are dealing 
with B9 and AO stars, and very few symmetric absorption lines 
are present in the spectra for the RV determination. Hydrogen 
lines are affected by emission, and several lines are affected by 
blending owing to rotation. For this reason, the independent es- 
timation of RV from the CO emission lines is also important. 

To apply the radial-velocity correction, we used the average 
value between of the three radial velocity determinations. We 
employed an RV of 16±2 km s' for HD 100546, and an RV of 
18±3kms-i for HD 97048. 

Absolute flux calibration was made by multiplying the tel- 
luric coiTected spectrum by the flux of a Kurucz model of the 
spectral type of the standard star at the observed wavelengths. 
The absolute flux calibration is accurate at the 20-30% level. 
Imperfections in the telluric correction and, most importantly, 
slit losses owing to the narrow slit and AO (Adaptive Optics) 
performance are the principal sources of uncertainty. The data of 
the settings Aref-2\\1 nm and /l,,,/=2123 run in the 2008 obser- 
vations were taken with the same exposure time (see Table A.l 1. 



Their combined spectrum was derived calculating the weighted 
average using the continuum flux as weight. The 2007 and 
2008 observations of HD 97048 have different exposure times 
(1920s and 320s respectively). The final combined spectrum of 
HD 97048 is the weighted average using the exposure time as 
weight. 

2.3. Position-velocity diagrams 

To construct position-velocity diagrams and to compare the 2D 
spectrum with 2D disk models, the 2D spectrum of the target was 
further processed: (/) it was corrected for the trace in the disper- 
sion direction by fitting a second-order polynomial to the PSF 
centers (obtained by a Gaussian fit) as a function of the wave- 
length, and then shifting each column of pixels such that the final 
trace is a straight line (correction on the order of 0.005 to 0.01 
pixel); (ii) to coiTect for the telluric absorption, each row in the 



^ Fiber-fed FEROS is the Extended Range Optical Spectrograph 
mounted at the ESO - Max Planck 2.2m telescope at la Silla, Chile. 
It covers the complete optical spectral region in one exposure. 
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spatial direction was divided by the ID spectrum of the standard 
star; (Hi) the 2D spectrum was re-centered in the wavelength di- 
rection such that the center of the line profile is at km s ' (in 
Sect. 3 we show that the center of the ID line profile is at a veloc- 
ity consistent with the rest velocity of the star. The re-centering 
correction of the 2D spectrum is smaller than the uncertainty of 
the radial velocity, which justifies this step); (iv) using the wave- 
length solution of the final extracted ID science spectrum, the 
2D spectrum was resampled in the dispersion direction such that 
wavelength scale had a uniform sampling; (v) the flux on each 
pixel was scaled in such a way that the extracted ID spectrum 
has the continuum equal to 1 ; finally, (vi) using the distance of 
the star and the CRIRES pixel scale of 0.086 arcsec/pixel, we 
derived the pixel scale in AU using as zero reference the pixel 
corresponding to the center of the PSF trace. In this way, a nor- 
malized, telluric corrected, wavelength and spatial calibrated 2D 
spectrum was obtained. 

We have two observations in two epochs with slits at differ- 
ent position angles for HD 97048. Because the line obsei-ved in 
2007 has a much better S/N than the line observed in 2008 (see 
Fig. l2|i, we employed only the data taken in 2007 to model the 
line profile and the analysis of the 2D spectrum (i.e., determina- 
tion of the spatial peak position (SPP), the FWHM of the PSF, 
and construction of the position-velocity diagrams). The 2007 
observations were taken with the slit in the N/S direction (see 



Table 2. 

limits. 



Summary of the measured line fluxes and 3cr upper 



Table A.l in the appendix). 

Two spectra were taken immediately after each other for 
HD 100546 : one with the wavelength setting centered at 2117 
nm and another centered at 2123 run. We produced an averaged 
2D spectrum by making a cut at -50 to +50 km s"' of the H2 1-0 
S(l) line in each observation and averaging the 2D frames using 
their continuum flux as weight. 

3. Results and analysis 

In Figs. [T] and |A.1| (in the appendix) we present the results of 
our survey. The H2 1-0 S(l) line is detected at the Scr level in 
HD 100546 and at the 8cr level in HD 97048. In the other 13 
targets the line is not detected. The H2 1-0 S(0) and 2-1 S(l) 
lines are undetected in all sources. We present in Table |2] the 
measured H2 line fluxes and upper limits. 

The H2 1-0 S(l) lines observed in HD 97048 and HD 100546 
are single-peaked. For HD 97048, a Gaussian fit to the line pro- 
file sets its center at -2.9+3 km s"' and the FWHM to 12.5 ± 1 
km s '. Because the line FWHM is broader than the spectral res- 
olution of CRIRES (~ 3 km s '), the line is spectrally resolved. 
In the 2007 2D spectrum, we observe that the FWHM of the ob- 
sei-ved PSF increases a few milli-arc-seconds at the position of 
the line, and we detect a clear motion of the spatial peak position 
(SPP, i.e., center of the gaussian fit to the observed PSF) along 
the line, which indicates that the emission is spatially resolved. 
The SPP displays the typical behavior for disk emission ruling 
out foreground emission (see Fig. [3]). 

For HD 100546, the Gaussian fit to the H2 1-0 S(l) line pro- 
file sets the center at 0.5+2 km s"' and its FWHM to 8 + 1 km 
s '. The line is spectrally resolved. In this source, the FWHM 
of the observed PSF slightly increases at the position of the line, 
which tentatively indicates that the emission is spatially resolved 
(see Fig. [3]). The SPP signal displays a shape suggesting an ex- 
tended emission component at positive velocities. However, be- 
cause the SPP signal is inside the 3cr noise level, it is difficult 
to establish at the S/N of our data whether this signal is real or 
noise residuals. The SPP signal observed is, nevertheless, consis- 
tent with the expected SPP signal of the rotating disk observed 





H2 1-0 S(l) 


H2 1-0 S(0) 


H2 2-1 S(l) 


Star 


^ =2121.83 nm 


A =2223.50 nm 


A = 2247.72 nm 


HD 58647 


<2.5 


<2.3 


<2.3 


HD 87643 


< 12 


<7.6 


<9.4 


HD 95881 


<3.5 


< 1.3 


< 1.1 


HD 97048 


9.6±2.9 


< 1.3 


< 1.5 


HD 100546 


5.4+1.6 


<2.1 


< 1.5 


HD 101412 


<4.6 


< 1.6 


< 1.7 


HD 135344B 


< 1.6 


< 1.7 


<3.5 


HD 141569 


<2.3 


<0.5 


<0.4 


HD 144432 


<0.2 


<0.2 


<0.3 


HD 150193 


<0.9 


<3.0 


< 1.8 


510ph 


<2.7 


<2.9 


<2.7 


HD 169142 


<0.9 


<0.5 


<0.7 


RCrA 


<7.2 


<9.4 


<6.6 


HD 179218 


<0.9 


< 1.0 


<0.9 


HD 190073 


<0.8 


<0.4 


<3.3 



Notes. All line fluxes and upper limits given in 10"''* erg s"' cm"'. The 
uncertainty on the integrated line fluxes corresponds to an uncertainty 
of 30% on the continuum flux. 



Table 3. Parameters of the 1-0 S(l) H2 fines detected. 



Star 



^0 
[nm] 



'•■center 

[nm] 



[km s 



FWHM 

'] [kms-'] 



HD 97048 
HD 100546 



2121.831 
2121.831 



2121. 81lH 
2121. 835H 



0.028 
0.014 



-2.9+3 
0.5+2 



12.5+1 
8+1 



Notes. *"' **' The error in the center and 6,1 is the combined error in the 
wavelength calibration of 0.5 km s ' and the uncertainty in the radial 
velocity of the sources. 



in CO (van der Plas et al. 2009) at the position angle of our ob- 
servations (PA~ 27°). Additional observations with higher S/N 
are required to test this. 

In Table [3] we summarize the parameters of the 1-0 S(l) H2 
lines detected in HD 97048 and HD 100546. Taking into account 
an uncertainty of 2-3 km s ' in the radial velocity of HD 100546 
and HD 97048, the emission is observed at a velocity consistent 
with the rest velocity of the stars. Because both sources have 
circumstellar disks, the most likely scenario is that the emission 
arises from the disk. The data could also be consistent with emis- 
sion of a disk-wind or outflow, but, in that case, it must have a 
low projected velocity (< 6 km s"^). 

To quantify the spatial extent of the emission lines detected 
in HD 97048 and HD 100546, we subtracted from the 2D spec- 
trum the average PSF observed in the continuum from -50 to 
-20 km s"' and +20 to +50 km s"' with respect to the line's 
position, and produced two diagrams. The first diagram is the 
position-velocity diagram of the 2D continuum-subtracted data 
(see Fig. [3]). In the second diagram, we summed all the counts 
from -15 to +15 km s"' and plotted the resulting cumulative 
counts as a function of the spatial off'set (see Fig.|4]i. 

In the case of HD 97048 these two diagrams show us that the 
H2 emission extends at least up to -150 AU (at negative veloci- 
ties) and +200 AU (at positive velocities). The position-velocity 
diagram displays the typical butterfly shape of a disk in rotation: 
the line emission at negative velocity tends to have a positive 
spatial offset, and the line emission at negative velocity tends to 
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Fig. 1. Obsei-ved CRIRES spectrum at the location of the H2 1-0 S(l), 1-0 S(0) and 2-1 S(l) Hnes in HD 97048 and HD 100546, 
the only two sources displaying H2 emission. Note that the H2 1-0 S(l) line in HD 97048 is the weighted average of the 2007 and 
2008 observations (see Fig.l2]). The spectra are presented at the rest velocity of the stars. 
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have a negative spatial offset. This is consistent with emission 
produced in a rotating disk and confirm the results suggested by 
the SPR 

The cumulative counts diagram of HD 100546 (see Fig. pi 
shows that the H2 1-0 S(l) line extends at least up to -50 AU to 
H-40 AU. The emission is slightly stronger at negative spatial off- 
sets. The position-velocity diagram displays the emission con- 
centrated at zero velocities and shows that the emission is much 
more compact (i.e., less extended) than in the case of HD 97048. 
At the S/N of the HD 100546 data, the signature of a rotating 
disk is not apparent, although weak extended emission is ob- 
sei-ved at positive velocities (see Fig. [3]). 

Note in Fig. |4] that for both stars the FWHM of the PSF at 
the position of the line (after subtracting the PSF of the con- 
tinuum) is larger than the FWHM of the PSF of the continuum 



(overplot in light gray in Fig.Bll. For HD 97048 the continuum- 
subtracted PSF FWHM at the H2 line position is ~520 mas, 
while the FWHM of PSF of the continuum is ~255 mas. For 
HD 100546 the continuum-subtracted PSF FWHM is -400 mas, 
while the FWHM of the continuum PSF is ~250 mas. This jus- 
tifies the approach of using Fig. HI to derive constraints in the 
spatial extension of the emission. Note that these are conserva- 
tive lower limits, because higher S/N data may reveal that the 
emission extends farther out. 

Both sources show evidence of extended envelopes from var- 
ious observations of gas and dust (e.g., Haitmann et al. 1993, 
Grady et al. 2001, in the case of HD 100546; Doering et al. 2007, 
Martin-Zaidi et al. 2010 in the case of HD 97048). However, the 
H2 emission observed is not consistent with emission from an 
envelope. In the case of HD 100546 the size of the envelope 
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Fig. 3. Observed and modeled spectrum for HD 97048 and HD 100546. For each star the upper three panels display the observed (in 
black) and the modeled (in red) H2 1-0 S(l) extracted ID line profile, the FWHM of the PSF, and the spatial peak position (SPP., i.e. 
center of Gaussian fit to the PSF). The lower three panels show the position velocity diagrams of the 2D spectrum after subtraction 
of the continuum PSF. The first panel displays the observed 2D spectrum, the second shows the modeled 2D spectrum, and the third 
their difference. We show a disk model with a =1.5, Rin-9 AU, and /?out=400 AU for HD 97048 and a disk model with or =1.0 
Rin^lQ AU, and /?out=80 AU for HD 100546. The spatial scale on the figures is 1 AU = 5.6 mas for HD 97048, and 1 AU = 9.7 mas 
for HD 100546. Note that the H2 lines have been shifted such that their center is at km s '. The 2007 and 2008 observations of 
HD 97048 have different slit PA. Because the 2007 data have much better S/N than the 2008 data, we considered for this figure and 
line modeling only the 2007 data. 



is ~1000 AU from dust scattering imaging (Grady et al. 2001). 
Because the line is spectrally resolved and its FWHM is much 
larger than the maximum expected Keplerian broadening at 1000 
AU (~ 3 km s"'), envelope emission is ruled out. In addition, if 
the emission would come from an envelope, its spatial exten- 



sion would be much larger than measured in the PSF FWHM. 
A similar set of arguments holds for HD 97048. Furthermore, in 
this source the shape of the SPP and the line position-velocity 
diagram rules out emission from an envelope, and an H2 line 
centered at the velocity of the star is inconsistent with the radial 
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represents the 2cr level of the background noise. The light gray 
line shows the scaled PSF of the continuum. 



velocity shift (14.5 km s"', Martin-Zai'di et al. 2010) obsei-ved 
in absorption lines (CH and CH^) tracing HD 97048's envelope. 

Finally, we note that the line flux measured in HD 97048 
is one order of magnitude higher than reported by Bary et al. 
(2008) (8.6+0.4 X 10"'^ erg s"' cm-^) using Phoenix (R ~ 6 
km s"') on Gemini South. The reason for this is unclear. It 
could be either due to observational reasons (difference in res- 
olution, i.e., with higher spectral resolution, it is easier to sepa- 
rate the line from the continuum; the use of AO also minimizes 
slit losses from the extended emission component) or the line 
might be intrinsically variable. Although our observations were 
not designed to test for variability, there is tentative evidence 
that the line is intrinsically variable. We obtained observations 
of HD 97048 with an identical telescope setup in two epochs 
separated by 15 months. In the first observation (April 2007) the 
line has a flux of 10 + 3 x 10"'"* erg s"' cm"^, and in the second 
epoch (July 2008) the line flux is 7 ± 2 x 10"'"* erg s"' cm"^ (see 
Fig. [2]). This change in the flux could be caused by variability 
of the line. However, note that the eiTor bars of the two mea- 
surements overlap. Therefore, the change seen in the line fluxes 
can also be due to the uncertainty on the determination of the 
continuum flux (30%). 

Comparing our measurements with those of Bary et al., we 
find that even allowing for a 30% uncertainty in the continuum 
fluxes, our line fluxes are still several times higher than those 
measured by Bary et al. (2008). This diff'erence cannot only be 
attributed to uncertainties on the continuum flux. Measurements 



especially designed for variability (e.g., spectroscopy and pho- 
tometry measurements simultaneously obtained) are needed to 
confirm the variability of the 1-0 S(l) H2 fine. 

3. 1. Line ratios and excitation mechanism. 

The line ratios of the H2 lines allow us to constrain the excita- 
tion mechanism and the temperature of the gas. Because the 1-0 
S(0) and 2-1 S(l) lines are not detected, we can only derive up- 
per limits to the H2 1-0 S(0)/l-0 S(l) and 2-1 S(l)/l-0 S(l) line 
ratios. For this the H2 1-0 S(0) and 2-1 S(l) line flux upper limits 
were divided by the measured H2 1-0 S(l) line flux minus 30% 
of uncertainty. In Table|4 we summarize our results. For compar- 
ison, we include in Table 4]the expected LTE emission line ratios 
for H2 at diverse temperatures, H2 emission excited by UV, X- 
ray, and shocks compiled by Mouri (1994) (see their Fig. 3 and 
Table 2), X-ray excitation models of H2 by Tine et al. (1997), 
and models of H2 emission from classical T Tauri star disks by 
Nomura et al. (2007) assuming dust grains with a spatially uni- 
form distribution, dust well mixed with the gas, and maximum 
grain radii Omax = lOfim. 

We find that in the case of HD 100546, the 1-0 S(0)/l-0 S(l) 
and the 2-1 S(l)/l-0 S(l) H2 line ratio upper limits exclude 
pure radiative UV fluorescent H2 emission from low-density gas 
(n < 10"* cm""*). However, the line ratio upper limits are too high 
to be able to distinguish between thermal, UV, X-ray, or shock 
emission. 

In the case of HD 97048, the 1-0 S(0)/l-0 S(l) and 2- 
1 S(l)/l-0 S(l) H2 line ratios exclude pure radiative UV flu- 
orescent H2 emission in low-density gas. Additionally, the 1- 
S(0)/l-0 S(l) line ratio also excludes emission by gas at 
T<2000 K, UV excited thermal and fluorescent H2 emission 
from dense gas (n > lO'* cm"^), and shock excitation mod- 
els. Measured H2 line ratio upper limits are consistent with the 
line ratios from thermal gas radiating at a temperature warmer 
than 2000 K and/or from H2 excited by X-rays. We note that 
HD 97048 was detected in X-rays by ROSAT by Zinnecker & 
Preibisch (1995). A short XMM-Newton observation (34 ksec; 
ObslD 0002740501, PI R. Neuhiiuser) serendipitously detected 
it again with a flux of about 8x10"'"^ erg s"' cm"^ (in the 0.3 - 8 
keV band; XMMXASSIST database, Ptak & Grifliths 2003). 



3.2. Line modeling witti Keplerian i\at-disk models 

We modeled the H2 line using a toy model that mimics emission 
from gas in a Keplerian orbit assuming a flat disk with known 
inclination and position angle (PA) (see Table B). The intensity 
of the emission is designed to decrease as I(R) oc (R/Ri^^y", with 
Rin being the inner radius, and R the radial distance from the star. 
We used a 2D grid from /?!„ to /?out and from to 360. For each 
grid point the emission is calculated by multiplying the intensity 
times, the solid angle times, a normalized Gaussian with FWHM 
of 3 km s"' (to simulate the spectral resolution) and center equal 
to the projected Keplerian velocity of the grid point. The disk 
emission is convolved in the spatial direction with a 2D Gaussian 
with the same FWHM as the observed AO point-spread function. 
Then, we overlaid a slit on the simulated image and produced a 
2D spectrum. The simulated 2D spectrum was rebinned in such a 
way that the pixel scale in the dispersion and the spatial direction 
is the same as in our CRIRES data. The rebinning was made con- 
serving the total amount of flux in the 2D spectrum. Finally, the 
flux of the simulated 2D spectrum was re-normalized such that 
the continuum has a value 1 in the ID extracted spectrum, and 
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Table 4. Measured upper limits to the H2 emission line ratios and 
theoretical line ratios expected for H2 at LTE, and H2 excited by 
UV, X-ray, and shocks (Mouri 1994, Tine et al. 1997, Nomura et 
al. 2007). 
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T = 2000 K, «H = 10^ - 10^ cm-3 
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X-ray irradiation 


0.23 


0.06 


UV irradiation 


0.25 


0.02 
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0.24 


0.03 



Notes. '"' Models of pure radiative UV fluorescent Hi emission spectra pro- 
duced in low-density (n < lO"* cm-^) cold isothermal photodissociation regions 
from Black & van Dishoeck (1987). '** Models of H2 infrared emission spectra 
in dense (n > lO"* cm"^^) static photodissociation regions exposed to UV radia- 
tion that both heats and excites the Ht gas from Sternberg & Dalgai'no (1989). 
(1) (b) Parajneters: nr = the total density of hydrogen atoms and molecules; x = 
UV-flux scaling relative to the interstellar radiation field; *'^' X-ray heating mod- 
els of Lepp & McCray (1983) assume an X-ray luminosity in the 1-10 keV band 
of 10^5 ergs s"'. Here are given the line ratios of their model b. *'" X-ray excita- 
tion models of Draine & Woods (1990) provide the expected H2 line emissivity 
efficiencies assuming a rate of absoiption of X-ray energy y = 2 X 1 0" " erg s" ' , 
"H = 10' cm"^, and monochromatic X-rays of 100 eV. In the case of the Ht 1-0 
S(l) line, they obtain efficiencies of 9.1 X 10^^7.9 X 10"^ and 4.8 X 10"' for 
an X-ray energy absorbed per H nucleus of 1, 10, and 62.3 eV respectively. In 
HD 97048 and HD 100546, we measured H2 1-0 S(l) luminosities of 1.9 X lO^** 
and 3.4 X 10^** erg s"' respectively. Using the Draine & Woods (1990) H2 1-0 
S( 1 ) line efficiencies, these H2 l-OS(l) luminosities would imply a Lx from 1 0'" 



n3l 



n29.5 , 



to 10 erg s" , somewhat brighter than the Lx ~ 10- measured in Herbig Ae 
stars (Telleschi et al. 2007). ''"' Tine et al. (1997) models do not prescribe a spe- 
cific X-ray input luminosity. They assume electrons of 30 eV and use as input 
parameters the ionization rate f (10^^ to 10"" s"'), nH(10 - lO' cm"'), and 
T = 500, 1000, 2000 K (see Tables 8 and 9 of Tine et al. 1997). These models 
include, in addition to the eft'ects of X-rays, the effects of collision processes of 
H2, with H2, H, and He on the resulting H2 emission spectrum. W^' Brand et al. 
(1989). <s' Nomura et al. (2007) CTTS models assume a Lx ~ lO'" ergs s"' 



and Lpuv ~ lO" ergs s '. 
maximum grain radii Umax = 



Here are given the line ratios of the models with 
lO^im. 



the amount of flux in the line inside ±20 km s and ±700 AU 
is equal to the flux of the line in the same region as the observed 
2D spectrum. This re-normalization was made to enable us to 
compare directly the normalized 2D observed spectrum with the 
2D spectrum computed with the model. 



Because the stellar parameters (mass, distance, inclination) 
are all constrained for HD 100546 and HD 97048 in the literature 
(see Table Island van der Plas 2009), the free parameters for the 
line modeling are a, Ri„, and /?ouf The exponent a parametrizes 
the dependence of the intensity as a function of the radius; it de- 
scribes the combined effect of the variation of the density, tem- 
perature, and emissivity of the disk. This assumption allows us 
to reduce the number of free parameters for the line modeling. 
For example, if we assume that the surface density is constant, 
an a = 2 would reflect an intensity directly proportional to the 
stellar radiation field. 

We compared the observed extracted ID line profile and the 
2D line spectrum (i.e., the 2D spectrum minus the continuum 
PSF) with their simulated counterparts for the parameter space 
0.1 AU < Rin < 450 AU, 10 AU < Rom < 500 AU, and < » < 
5.5, calculating for each model the reduced ;^f^ statistic; 
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Ared 
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Here, O are the observed data points and M the modeled values. 
The variance cr^ was determined from the continuum emission 
between -50 and -20 km s"', and +20 and +50 km s"' at -300 
to - 100 AU and +100 to +300 AU. x^s-a ^^^ calculated over the 
data between -20 and +20 km s"'. A^ is the number of valid data 
points between -20 and +20 km s"'. The number of degrees of 
freedom (A^ - 4) is obtained because we use three free parame- 
ters for the fit (/?in, Roui, a), and one degree of freedom is taken 
because the simulated line is re-normalized to match the total 2D 
flux of the observed line. 

3.2.1. Modeling results 

In Fig. pi we show x^^.^ contours for the continuum-subtracted 
2D spectrum for diverse values of Rom in the a vs Ri„ parame- 
ter space, and the;^'^^^ contours for the ID extracted line profile 
as a function of Ria for different values of Rom for a fixed a. 
In the following paragraphs we describe the modeling results in 
each source individually. In TablelSlwe present a summary of the 
modeling of our data. 

HD 97048: The 2D spectrum clearly shows an extended 
emission component extending at least up to 200 AU. The x^^^ 
contours for models with Rom > 200 AU indicate that the solu- 
tions with minimum ;^'j^g^ converge to R^^ < 10 AU and a varying 
1.0 to 1.5 (see Fig.pb. The;i'^g^ decreases by increasing Rom- This 
suggests that the emission may extend at least to 300-400 AU. 
To constrain the multiple values of Ri^ giving similar x^^.,}^ in the 
2D spectra, we used the fit to the extracted ID line profile. By 
fixing a to 1.5, we calculated ;t^j^^j for the line profile for diverse 
Rom > 200 AU. Because the line observed is single peaked, so- 
lutions with radii larger than 3 AU provide a better fit. The x^^^ 
of the line profile decreases up to Ri„ ~ 8 - 9 AU, then the x^^^ 
starts to rapidly increase (see Fig.l3]top panel, upper right). This 
suggests that the Rin that best describes the line shape and the 2D 
spectrum simultaneously is ~ 8 AU. In summary, flat disk mod- 
els indicate that the H2 emission observed in HD 97048 is most 
likely produced in the disk starting at 5 < Rin <10 AU, and that 
it extends to a Rom of a few hundred AU. In Fig. |3] we present 
the line profile, FWHM, SPP, 2D spectrum, 2D model spectrum 
and their difference for a disk model with /?in=9 AU, Rom-'^OQ 
AU and «= 1.5. 

HD 100546: The detected H2 1-0 S(l) in HD 100546 is weak 
and the amount of counts in the 2D spectrum is low. Therefore, 
the direct comparison to 2D model data is challenging. The first 
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Fig. 5. Reduced ;i'^ contours of the H2 2D and ID line profile modeling for HD 97048 and HD 100546. For each star theirs? three 
panels display the reduced x^ contours for the 2D spectrum as a function of a and /?!„ for diverse values of /?out ■ The value of the 
minimum ;fj^j is indicated in the title. The contour in solid black is the contour of the minimum ;i^j^^^ plus 0.1. Subsequent contours 
in gray dotted lines are at the minimum ;t'^g^ plus 0.5, 1, 3, 5, 10 in the case of HD 97048, and at the minimum plus 0.2 and 0.5 in 

the case of HD 100546. The rightmost panel shows the ID line profile ;^'^g^ as a function of /?!„ for several values of 7?out for a fixed 
a. For HD 97048 a'=1.5 is used, in the case of HD 100546 ^=1.0 is employed. 



observational fact is that the line observed is single-peaked and 
that the line wings are < 10 km s"'. This is better described 
by disk models with emission at large inner radii /?in 30 -100 
AU (see lower right panel of Fig. B}. In contrast, the position 
velocity diagram of the 2D spectrum displays the emission rela- 
tively concentrated at /?out <50 AU, which favors disk solutions 
at 7?in < 20 AU (see lower left panels of Fig. l5]). Thus, the best 
models are a compromise between reproducing the single peak 
and the lack of high velocity wings in the line profile (i.e., favor- 
ing a large R^), while keeping the emission not too extended to 
fit the 2D spectrum. In the lower-left panels of Fig. [5] we present 
the;if^g^ contours for flat disk models with 7?out - 50, 80 and 100 
AU. We observe that there is a large family of solutions with 
similar ;t'^g^. If we assume a to be larger than 1.5, the solutions 
converge to a disk with R\n between 5 and 20 AU. If a is lower 
than 1.5, we can only say that R\^ should be smaller than 20. To 
break this degeneracy, the ID line profile is useful. In the lower- 
right panel of Fig. pi we display x^^^^ of the ID line profile as a 
function of R\^ for aiverse values of /?out- This plot shows that 
the x^^^^ is almost constant up to 10 AU, then decreases reach- 
ing its minimum between 30-100 AU (depending of /?out), and 
then increases rapidly. This plot suggests that /?[„ > 10 AU, and 
that for a given a the best combined fit of line profile and 2D 



spectrum is given by the maximum /?;„ inside the contour of the 
minimum ;(f^gj in the 2D spectrum. In summary, flat disk models 
indicate that the H2 emission observed in HD 100546 is most 
likely produced in the disk starting at 10 < 7?in < 20 AU and 
that the emission extends at least to 50 AU. In Fig. [3] we present 
the line profile, FWHM, SPP, 2D spectrum, 2D model spectrum, 
and their difference for a disk model with R\^—2^ AU, i?out=80 
AUanda^LO. 

The flat disk-models used here are "toy models" and are 
oversimplifications of the real structure of the disks, which are 
indeed flared-disks with (at least in the case of HD 100546) a 
complex inner disk structure in the dust (see Section 4.2). In a 
flared-disk geometry the outer regions of the disk are more ex- 
posed to radiation than in a flat-disk geometry. Therefore, it is 
expected that the effect of the flared geometry would be to dis- 
place the x^^^^ contours towards lower values of alpha to account 
for the increment in the flux contributions from the outer disk. 
Since the inner radius is less affected by the flaring, the inner ra- 
dius of the emission of a flared disk will not change considerably 
with respect to the inner radius of a flat disk. Note in particular 
that the lowest values f^^ x\^a arc >1.5, which formally means 
that the models do not fit the data correctly. Nevertheless, the 
main message that the flat disk models are telling us is that most 
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of the H2 emission observed is produced at R »5 AU. This 
is an important message because, in the context of passive disk 
models with gas and dust in thermal equilibrium developed for 
Herbig Ae/Be stars (e.g., Dullemond et al. 2001), the H2 gas at 
temperatures of few thousand K, which is necessary for emitting 
the near-IR lines, is expected to be located in the inner parts of 
the disk only up to a few AU. We will discuss this in detail in 
Section 4. 



4. Discussion 

4.1. H2, CO, and [Oi] emitting regions. 

As a summary of our observations and modeling, we can con- 
clude that in both sources the observed near-IR H2 emission ap- 
pears to be produced at distances larger than 5 AU. In the case 
of HD 97048 the emission extends at least up to 200 AU and in 
the case of HD 100546 the emission extends at least up to 50 
AU. In both targets, two other gas tracers, the [Oi] line at 6300 
A and the CO ro-vibrational emission band at 4.7 /zm, have been 
detected (e.g. van der Plas et al. 2009, Acke et al. 2005). In addi- 
tion, in the case of HD 97048 the H2 0-0 S(l) Une at 17 fim has 
been observed (Martin-Zaidi et al. 2007). How does the near-IR 
H2 emitting region compare with these other tracers? 

Seeing-limited observations of the 0-0 S(l) H2 line at 17 yum 
by Martin-Zaidi et al. (2007) set an upper limit to the extension 
of the emission of 35 AU. Because the 17 yum line is spectrally 
unresolved in the Martin-Zaidi et al. spectra (FWHM~30 km 
s"'), the profile cannot be directly compared with our CRIRES 
1-0 S(l) profiles. It can only be said that the inner radius pro- 
ducing the 17 jum line should be larger than a few AU, otherwise 
the line would have been spectrally resolved. Given that the 1-0 
S(l) H2 emission at 2 /im is located farther out in the disk (at 
least up to 200 AU), the two data sets clearly indicate that the 
emitting regions of the 2 and 17 /im line are different, with the 
1-0 S(l) Hne excited farther out than the 0-0 S(l) line. Note that 
in the case of thermal disk emission, the contrary is expected, 
namely that the 17 fim line is produced farther out because the 
temperature decreases with the radius. This indicates that the ex- 
citation mechanism of the two lines is different, or at least very 
inefficient at producing detectable levels of 0-0 S(l) at large dis- 
tances. Models of H2 emission from disks around T Tauri stars 
(e.g. Nomura et al. 2007), predicted 0-0 S(l)/ 1-0 S(l) line ratios 
raging from 0.13 to 40 for fixed grain sizes (lOjum to 10 cm) and 
30 to 70 for models including grain coagulation and settling. The 
measured line ratio from our and Martin-Zaidi's observations is 
0.25^gj assuming a 30% uncertainty in the continuum. This ra- 
tio is similar to that suggested for the UV and/or X-ray excitation 
models with grains of 10 fim size. It does not agree with the ra- 
tio from models with larger grains or those implementing dust 
coagulation and settling. 

In the case of the [Oi] line at 6300 A and the CO ro- 
vibrational spectra at 4.7 fj.m, the spectral resolution of the data 
is similar to that of our observations. Therefore, the line profiles 
can be used to compare the different emitting regions. In Fig.l6] 
we plot the [O i] line spectrum, the composite y = 1 - CO ro- 
vibrational spectrum, and our H2 1-0 S(l) line observations. To 
facilitate the comparison, the data are presented in such way that 
the lines are centered at the rest velocity, the continuum is set to 
zero, and the peak of the lines is set to 1 . All spectra have simi- 
lar spectral resolution. In this way we can compare the width and 
the high-velocity wings of the different transitions. Additionally, 
we include in Table |5] the results of kinematic modeling of the 



Table 5. Modeling results summary. 



Star 


Species 


Model 


[AU] 


"out 

[AU] 


a 


HD 97048 

/ = 47° 
PA=175° " 


H2 

CO" 

[Oi]" 


Flat disk 
Flat disk 
Flared 


5-10 
11 
0.8 


>200 
>50 
-50 


1.5 
2.5 


HD 100546 

i = 43° 
PA=160° 


H2 

CO" 

[Oi]" 


Flat disk 
Flat disk 
Flared 


10-20 

8 

0.8 


>50 
>50 
-50 


ind. 

2.5 



Notes. <"' van der Plas (2009); <*' Acke et al. (2005); <''' The posi- 
tion angle (PA) of the major axis of the disk around HD 97048 on the 
sky is 175±1° E of N, derived from the 8.6yum image of Lagage et 
al. (2006). This value agrees with the PA previously determined from 
spectro-astrometry of the [O i] 6300 A line (160+19°, Acke et al. 2005). 



u = 1 - CO (van der Plas et al. 2009) and [OI] (Acke et al. 
2005) line profiles. 

As a general trend in both sources, we observe the follow- 
ing characteristics: (i) The [Oi] line is the broadest line and has 
velocity wings extending to ±50 km s"'; (ii) The CO line pro- 
file is narrower than the [Oi] line and broader than the H2 line; 
(Hi) The [Oi] profiles are broad and double-peaked; (iv) The 
CO lines are broad, and the case of HD 97048 they are consis- 
tent with a double-peaked profile (van der Plas (2010) detected 
the L- = 2 - 1,3 - 2, and 4 - 3 CO transitions in HD 100546 
and HD 97048, the observed line profiles are consistent with 
double-peaked profiles). In contrast, the H2 line is single-peaked. 
Assuming that we are observing disk emission, these character- 
istics indicate that each gas tracer is produced in a different radial 
region of the disk. The [O i] line is produced the closest to the 
star, it is followed by the CO emission at a larger distance, and 
finally the H2 line farther out in the disk. 

4.2. H2 near-IR emission and observations of dust in the 
inner disk 

HD 97048 and HD 100546 have been extensively studied in 
dust continuum emission, including studies with ground-based 
infrared interferometers in the case of HD 100546 (e.g., Benisty 
et al. 2010). How does the H2 near-IR line-emitting region com- 
pare with inner disk structure deduced from dust observations? 

In the case of HD 100546, recent simultaneous modeling of 
the spectral energy distribution (SED) and near-IR interferome- 
try (Benisty 2010, TatuUi et al. 201 1) suggested a dust disk struc- 
ture consisting of a tenuous inner disk located around 0.25 AU 
to 4 AU, followed by a gap devoid of dust, and a massive outer 
disk starting at 9 to 17 AU that extends up to a few hundred AU. 
Comparing this disk dust structure and the H2 near-IR emitting 
region, we find that the inner radius of the region emitting the H2 
line and the inner radius of the massive outer dusty disk are re- 
markably similar. If the inner radius of H2 near-IR emission and 
the beginning of the outer thick disk are indeed the same, this 
would suggest that the presence of the gap in the dust is favoring 
the excitation of H2 at large distances from the star. However, it 
still remains to be explained in this scenario, why [OI] and CO 
emission are present at radii smaller than H2. 

In the case of HD 97048, no near-IR interferometry obser- 
vations are reported in the literature. The most comprehensive 
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Fig. 6. Observed line profiles of H2 1-0 S(l) emission (in gray, 
histogram), combined y = 1 -0 CO emission, and [Oi] emission 
for HD 100546 (top) and HD 97048 (bottom). All three spectra 
have similar spectral resolution. 



study of the dusty disk structure is from Doucet et al. (2007), 
who simultaneously modeled the SED, mid-IR spectra and spa- 
tially resolved mid-IR imaging in the continuum and PAH fea- 
ture at 11.3 i-im. Those authors suggest a continuos flared disk 
model with a puffed-up inner rim at radius 0.4 AU, outer radius 
370 AU, and flaring index 1.26. The extension of the dust emis- 
sion down to a fraction of AU combined with the detection of gas 
tracers such as the [OI] line in HD 97048 also down to a fraction 
of AU (Acke et al. 2005), contrasts with the H2 near-IR emitting 
region at a radius larger than 5 AU found in our observations. 
The reason for this discrepancy is unclear. 

4.3. H2 near-IR emission in the context of diverse disk 
structure models 

In summai-y, we find fliat in HD 100546 and HD 97048, H2 near- 
IR emission appears only to be present at large radii. But, as 
discussed, within these radii material is still present within the 
disk, as shown by the presence of warm dust (e.g. Bouwman et 
al. 2003; TatulH 2011; Doucet 2007), the kinematics of the [Oi] 
emission line at 6300 A (Acke & van den Ancker 2006) and the 
CO ro-vibrational band at 4.7 yum (van der Plas et al. 2009). The 
situation as seen in the H2 1-0 S(l) line reported here is some- 
how analogous to what is seen in CO around these and other 
Herbig stars (Brittain et al. 2009; van der Plas et al. 2009; 2010), 



where there also appears to be a deficit of CO emission com- 
ing from the inner regions of the disk. Van der Plas et al. (2009) 
have argued that this may imply that CO is efliciently depleted 
close to the star. However, it is unlikely that this explanation will 
also hold for H2, given its much greater abundance; one would 
only expect H2 to be depleted in a disk that is gas poor overall. 
Moreover, the detection of [O i] emission close to the star sug- 
gests that the disks around HD 97048 and HD 100546 are still 
gas-rich, which eliminates this possibility. To address this and 
explain the numerous non-detections of H2 near-IR emission, 
we discuss below our results in the frame of passive disks mod- 
els with thermally coupled gas and dust, passive disks models 
with thermally decoupled gas and dust, and disk photoevapora- 
tion. 



4.3.1 . H2 near-IR emission in tine context of passive disks 
witli tliermally coupled gas and dust 

Standard passive disks models of Herbig Ae/Be disks, i.e., mod- 
els in which gas and dust are thermally coupled in the disk sur- 
face layer (e.g., Dullemond et al. 2001), show that in the outer 
regions of the disk {R > 5 AU) the temperatures in the disk are 
too low to emit efficiently in the detected H2 1-0 S(l) line, which 
requires gas of a few thousand degrees to be excited. Employing 
Chiang and Goldreich (1997) Herbig Ae two-layer disk models 
(using the implementation by Dullemond et al. 2001 described 
in Carmona et al. 2008b), we calculated the amount of gas in the 
surface layer at Tsmf >1000 K for disks of mass 10, 40 and 100 
Mj (see Figs. 4 and 5 of Carmona et al. 2008b). We found that 
the mass of gas Tsurf >1000 K is ~0.01 M,„oon- Therefore, in the 
inner disk, where temperatures are sufficiently high to thermally 
(i.e., by collisions with dust) excite the H2 lines, the amount of 
gas in the optically thin surface layer of the disk is small, quench- 
ing the formation of the near-infrared H2 lines. Because CRIRES 
observations are typically sensitive from 0.1 to 1 Mmoon of H2 
gas at T>1000 K (Carmona et al. 2007), our first conclusion is 
that the lack of near-IR H2 emission in the 13 sources with non- 
detections fully agrees with what would be expected from pas- 
sive Herbig Ae/Be disks with thermally coupled gas and dust. 

The puzzling new observational result presented in this pa- 
per is that H2 emission is present at radii >>5 AU in the disks 
around the Herbig Ae/Be stars HD 97048 and HD 100546 
Because this is not expected in the context of passive disks with 
thermally coupled gas and dust, we conclude that at least in these 
two objects, the disk structure is not accurately described by sim- 
ple passive disk models. Either (/) at large radii molecular gas 
and dust in the disk surface layer have departed from thermal 
coupling (i.e., Tgas > Tjust) and/or (//) the disk atmosphere must 
be much more extended at large radii, for example in the form a 
disk wind caused by photoevaporation. We note that there is ob- 
servational evidence that dust and gas thermally decouple in the 
surface layers of Herbig Ae/Be disks. Fedele et al. (2008) stud- 
ied the dust (traced by the 10 fiva feature) and the gas (traced by 
the [O i] line at 6300 A) in the disks of a sample of Herbig Ae/Be 
stars. They found a difference in the gas and dust vertical struc- 
ture beyond 2 AU in the Herbig Ae/Be star HD 101412, thereby 
providing evidence of gas-dust decoupling in a protoplanetary 
disk atmosphere. 



'^ Detections of near-IR H2 emission displaying an extended ring 
structure at 73-219 AU away from the star have been reported toward 
the weak-lined T Tauri star (WTTS) DoAr 21 (O. Panic 2009 PhD the- 
sis). DoAr 21 is one of the brightest X-ray WTTS (Neuhauser et al. 
1994). 
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4.3.2. H2 near-IR emission in tine context of passive disks 
witli tliermally decoupled gas and dust 

Models of Herbig Ae/Be hydrostatic disks allowing the decou- 
pling of gas and dust in the disk surface layers and including 
additional heating mechanisms of the gas have been recently 
developed (e.g., ProDiMo, Woitke et al. 2009). These types of 
models can provide a hint at the possible origin of near-IR H2 
emission observed in HD 100546 and HD 97048. In these mod- 
els, the temperature of the gas in the surface layer of the disk 
can be much higher than the gas temperature in the passive disk 
models where gas and dust are thermally coupled. The thermally 
decoupled gas reaches temperatures up to a few thousand K at 
distances up to 10-20 AU in Herbig Ae/Be flared disks. In self- 
shadowed disks, the temperature is also high, but is constrained 
to a radial extent of < 5 AU (van der Plas, 2010). 

In the ProDiMo model, the transition between H and H2 in 
the disk is determined by the balance between the formation rate 
of H2 on grains and the dissociation of H2 by the stellar UV 
field. At high temperatures H2 has a low formation efficiency. 
The sticking of H on grains to form H2 is high, but the lifetime 
of H on the surfaces is very short owing to thermal desoiption 
if the grain temperature is also high. As soon as the H2 forma- 
tion rate can balance the H2 dissociation, abundances start to 
build quickly. Otherwise abundances stay low (on the order of 
10"^). In the models discussed by van der Plas (2010), the H/H2 
transition occurs at the location where the gas temperature drops 
below < 1000 K, and/or where the UV opacity is high enough to 
protect the H2 against photo-dissociation, or around 10 AU. 

Thermally decoupled gas and dust passive disk models 
would naturally explain why the H2 emission is observed in flar- 
ing disks but not in self-shadowed disks, and explain the numer- 
ous non-detections. However, it is not clear why then H2 near- 
IR emission is not present in all the sources with flared disks 
in our sample. One possible explanation is that HD 97048 and 
HD 100546 are flared disks with the earliest spectral types (AOV 
and B9V respectively) of the flaring disks that we have observed. 
The earlier the spectral type, the more stellar radiation is avail- 
able to heat the upper-disk layers. An additional challenge for 
the hydrostatic thermally decoupled gas-dust models is that at 
least in the case of HD 97048, the line ratio limit of the detected 
1-0 S(l) line to the (undetected) 1-0 S(0) line sets the gas tem- 
perature around 2000 K. Because gas at T ~ 1000 K is more 
abundant than gas at T ~ 2000 K, copious amounts of H2 emis- 
sion from gas atT ~ 1000 K would be expected, and gas at this 
temperature is not observed. This suggests that in addition to de- 
coupling of the dust and gas some change in the disk structure at 
radii >5 AU is also necessary to explain our observations. 

4.3.3. An extended disk atmosphere caused by 
photoevaporation ? 

It is interesting to consider the effect of photoevaporation in the 
context of the possible presence of an extended disk atmosphere 
at large radii in the cases of HD 97048 and HD 100546. H2 near- 
IR emission is detected at a velocity consistent with the rest ve- 
locity of the stars. However, owing to uncertainties on the radial 
velocity of the sources, the possibility that the lines have a blue 
shift of a few km s"' cannot be completely excluded. A small 
blue shift would be consistent with an origin of the emission in 
a photoevaporative disk wind. 

Photoevaporation of disks by ionizing photons is a conceptu- 
ally simple process in which energetic radiation (far-UV (FUV), 
6 eV < /iv < 13.6 eV; extreme-UV (EUV), 13.6 eV < hv < 0.1 



keV ; and X-rays hv > 0. 1 keV) from the central star heats hy- 
drogen at the disk surface, producing a hot ionized layer. Near 
the star this ionized zone above the disk is expected to be almost 
static. However, at large radii from the star the thermal layer will 
be unbound, powering a theiTnally driven disk wind (Shu et al. 
1993; Hollenbach et al. 1994). The critical radius at which the 
upper layer of the disk will become gravitationally unbound is 
on the order of 10 AU in the case of T Tauri stars and scales 
linearly with the mass of the central star (Font et al. 2004). This 
critical radius is compatible with the inner radii for the H2 1-0 
S(l) emission observed in HD 97048 and HD 100546. 

The H2 1-0 S(l) emission observed in both HD 97048 and 
HD 100546 could thus be due to the presence of an extended 
disk atmosphere (i.e., disk wind) at large radii owing to photoe- 
vaporation. 

In the case of HD 97048, this low velocity disk wind could 
also be responsible for the pure rotational H2 emission at 17 //m 
that has been previously detected (Martin-Zaidi et al. 2007). A 
critical radius larger than 10 AU could explain why the 17 jum 
line is observed spectrally unresolved. However, it remains to be 
understood, why the emission at 17 //m appears to be confined 
at R< 35 AU, while the emission at 2 /jm micron extends up to 
hundreds of AU. One possibility is variability of the 17 //m and 
2 /urn line because the two data sets were not obtained simulta- 
neously. 

In the case of HD 100546, Lecavelier des Etangs et al. (2003) 
and Martin-Zaidi et al. (2008a) reported H2 in absorption in the 
FUV domain with Ff/5/iPlat a radial velocity consistent with a 
circumstellar origin at the resolution of FUSE (R~15000). Warm 
and hot gas at 760 K and 1500 K kinetic temperature were 
measured and revealed excitation conditions of H2 clearly dif- 
ferent from those observed in the interstellar medium. Because 
the FUSE line of sight does not pass through HD 100546 disk, 
those authors concluded that the H2 observed by FUSE is not 
located in the disk and suggested that the H2 absorption might 
be produced by an FUV-driven photoevaporative wind from the 
outer parts of the disk. An interesting possibility is that precisely 
this hot H2 at T>1000 K from a photoevaporative wind seen by 
FUSE is being traced by the near-IR H2 emission lines. 

Models of EUV photoevaporation of disks around T Tauri 
stars predict for disks with inclination ~45° single-peaked emis- 
sion lines of ionized gas with blueshifts ranging from -12 to -7 
km s ' and line FWHM ranging from 26 to 30 km s ' (e.g., 
[N 11], [S II] Font et al. 2004; [Ne 11] Alexander et al. 2006, 2008). 
These values are somewhat higher than the line center and line 
FWHM of the H2 lines observed in HD 100546 and HD 97048 
(see Table l3]l, even taking into account the uncertainty on the 
radial velocity. If the H2 emission observed is linked to an EUV 
photoevaporative wind, the velocity of the outflowing gas should 
be lower than a few km s"' to be consistent with our observa- 
tions. The smaller line FWHM observed may be explained by 
the fact that the H2 line is dominated by contributions at large 
radii. 

In the context of photoevaporative disk winds driven by FUV 
or X-rays (e.g., Gorti & Hollenbach 2009, Ercolano & Owen 
2010), there might be a connection between the presence of ex- 
tended H2 near-IR emission and the rich and spatially extended 
PAH spectrum observed in HD 97048 and HD 100546 (van 
Boekel et al. 2004, Habart et al. 2006). In contrast to EUV-driven 
disk winds, FUV and X-ray driven disk winds are predominantly 
denser, cooler, mainly neutral, and may drag along small grains. 



Far Ultraviolet Spectroscopic Explorer. 
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The photoevaporative wind scenario does immediately raise 
the question why the H2 lines we observe here have only been 
detected in HD 97048 and HD 100546 and not in the other 13 
Herbig Ae/Be stars in which we searched for the emission. In the 
case of objects with self-shadowed disks, the temperatures in the 
disk surface layer may not be high enough to launch a disk wind. 
In the case of flared disks, HD 97048 and HD 100546 are of rel- 
atively early spectral types (AO and B9, respectively), so pho- 
toevaporation or direct disk heating from stellar UV radiation 
could be more prominent in these stars. In addition, HD 97048 
and HD 100546 are the Herbig stars with the most extreme UV 
fluorescence detected in CO ro-vibrational emission at 4.7 yum 
(see van der Plas et al. 2009, 2010) and have the strongest PAH 
features at 8.6 and 1 1 yum (Acke & van den Ancker 2004). The 
two other sources with flared disks in our sample, HD 169142 
and HD 179218, have slightly different physical properties than 
HD 97048 and HD 100546. HD 169142 has a spectral type A5. 
Because it is fainter than HD 97048 and HD 100546, the gas 
in the disk may be not sufficiently excited to produce detectable 
levels of H2 near-IR emission. HD 179218 has a much higher lu- 
minosity than the other Herbig stars in the sample, but happens 
to have similar Teff as HD 100546. The higher luminosity might 
cause photodissociation of H2 to larger distances from the star. 
In this object van der Plas (2010) have detected warm CO at 4.7 
yum likely excited trough UV fluorescence, albeit much weaker 
than in HD 97048 and HD 100546. The fact that HD 179218 
has a distance of 240 pc may explain why the CO emission is 
weaker and the H2 emission undetected. Alternatively, the pho- 
tons responsible of photoevaporation could be driven by accre- 
tion, which may be somewhat above average in HD 97048 and 
HD 100546 (van den Ancker 2005). 

5. Conclusions 

We presented CRIRES high-spectral resolution observations of 
H2 ro-vibrational emission at 2 /im toward a sample of 15 Herbig 
Ae/Be stars. We detected 1-0 S(l) H2 emission from HD 97048 
and HD 100546 at a velocity consistent with the rest velocity 
of the stars; in the other sample stars the emission is absent. 
None of the sources displays the H2 1-0 S(0) or the 2-1 S(l) 
line. In the case of HD 97048 the emission extends at least up 
to 200 AU. The spatial peak position (SPP) and the position- 
velocity diagram show that the emission at negative velocities 
has a positive spatial offset and that the emission at positive 
velocities has a negative spatial offset. This provides evidence 
that the emission is observed from rotating material, most likely 
linked to the circumstellar disk. The H2 1-0 S(l) flux measured 
in HD 97048 is one order of magnitude stronger than the line 
flux reported by Bary et al. (2008) in observations carried out 
in 2003. Because the difference in the line fluxes between our 
and Bary et al.'s (2008) measurements is much larger than the 
uncertainty on the continuum flux, this suggests line variability. 
In the case of HD 100546, the emission extends at least up to 
50 AU. The SPP is compatible with the signature of material in 
rotation; however, observations with higher S/N are required to 
firmly establish this. 

For HD 97048 and HD 100546, we determined upper limits 
to the 1-0 S(0)/l-0 S(l) and 2-1 S(l)/l-0 S(l) line ratios and con- 
strained the excitation mechanisms. In both sources pure radia- 
tive UV fluorescent H2 emission from low-density gas (n < 10"* 
cm"^) is excluded. In the case of HD 100546 the line ratios upper 
limits did not permit us to distinguish between UV (thermal and 
fluorescent) excited H2 emission in dense gas (n > 10^ cm"^). 
X-ray, shocks, or thermally excited H2 emission. In the case of 



HD 97048, the line ratios upper limits are consistent with ther- 
mal emission at T >2000 K and/or H2 gas excited by X-rays. 

Assuming a flat disk in Keplerian rotation, we modeled the 
line profiles, the PSF FWHM, the SPP and the 2D spectrum of 
the 1-0 S(l) H2 fine observed in HD 97048 and HD 100546. In 
both cases the models indicate that the emission arises at large 
(»5 AU) radii. In the case of HD 97048, the models suggest 
that the emission is produced starting from a /?[„ of 5-10 AU and 
extends at least to 200-300 AU. Solutions converge to an a ~1.5. 
In the case of HD 100546, the models suggest that the emission 
is produced starting from a R^^ of 10-20 AU and extends at least 
to 50 AU. 

In conclusion: 

a) The non-detections of near-IR H2 can be explained (/) as 
a natural consequence of a passive disk structure in which the 
gas and dust are thermally coupled (i.e., the amount of optically 
thin hot H2 gas in the inner disk is too low to produce detectable 
levels of H2 emission). Or (//) in the context of a passive disk 
with thermally decoupled dust and gas; here the H2 close to the 
star is photo-dissociated, and temperatures in the outer disks of 
self-shadowed disks are too low to produce H2 emission. 

b) Most of the H2 emission observed in HD 100546 and 
HD 97048 is produced aX R » 5 AU and extends at least to 
R > 5Q - 200 AU. This location is much farther out in the disk 
than the radii up to a few AU where temperatures in the disk 
are expected to reach the 1000-3000 K that is needed to excite 
the observed H2 lines thermally in standard passive disk mod- 
els with thermally coupled gas and dust. The emission can be 
explained by (/) disk models with thermally decoupled gas and 
dust that allow for high gas temperatures at large radii in flaring 
disks. At small radii, H2 is photo-dissociated by the strong stellar 
UV field and H2 formation efficiencies are low. After ~ 10 AU, 
H2 formation becomes more efficient at balancing photo dissoci- 
ation, thus leading to detectable H2 gas in the disk atmosphere; 
or (//) by an extended disk atmosphere owing to a photoevapo- 
rative disk wind. More detailed modeling would be of great help 
to understand the origin of the H2 near-IR emission observed in 
HD 100546 and HD 97048. 
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Table A.l. Summary of the observations 



[nm] 



2117.7 
2123.4 
2236.1 

2117.7 
2123.4 
2236.1 



Date UT start 

dd/mm/yyyy hh:mm 



'■exp 
[S] 



slit PA 

n 



PSF FHWM 

SCIENCE STD 

[mas] [mas] 



HD 58647 (calibrator: HD 44402) 

2117.6 27/02/2008 03:23 160 1.1 118 220±3 207±5 

2123.3 27/02/2008 03:29 160 1.1 117 219±3 218±5 

2236.1 27/02/2008 03:34 160 1.1 116 209+3 207±5 

HD 87643 (calibrator: HD 81188) 

2117.6 14/03/2008 04:07 32 1.2 20 232±5 497±24 

2123.3 14/03/2008 04:12 32 1.2 22 225±5 387±36 
2236.1 14/03/2008 04:17 32 1.2 24 223±8 721±24 

HD 95881 (calibrator: HIP 061 199) 

2117.7 02/07/2008 23:06 160 1.6 39 236±4 256±8 

2123.4 02/07/2008 23:11 160 1.6 40 214±5 278±9 
2236.1 02/07/2008 23:53 160 1.7 52 274±45* 237+3 



HD 97048 (calibrator: HIP 061 199) 
2117.7 04/07/2008 23:20 160 1.8 40 250±5 

2123.3" 01/04/2007 03:35 1920 1.7 278±4 



256±8 



2123.4 


04/07/2008 23:25 160 1.8 41 


226±8 


278+9 


2236.1 


05/07/2008 00:17 240 2.0 122 


237+41 


237+3 




HD 100546 (calibrator: HIP 061585) 






2117.7 


03/07/2008 22:50 160 1.5 26 


253+3 


209+3 


2123.4 


03/07/2008 22:55 160 1.5 28 


257+3 


202+3 


2236.1 


03/07/2008 23:40 160 1.5 41 


224+6 


228+3 




HD 101412 (calibrator: HD 106490) 






2117.6 


15/03/2008 06:42 240 1.3 41 


214+6 


540+40 


2123.3 


15/03/2008 06:47 240 1.3 43 


206+7 


620+31 


2236.1 


15/03/2008 06:55 240 1.3 46 


219+30 


441+16 



HD 135344B (calibrator: HD 143118) 

2117.6 23/02/2008 07:44 240 1.1 -73 

2123.3 23/02/2008 07:50 240 1.1 -72 
2236.1 23/02/2008 07:56 240 1.1 -70 

HD 141569 (calibrator: HIP 061 199) 

2117.7 03/07/2008 02:47 160 1.1 146 

2123.4 03/07/2008 02:54 160 1.1 143 
2236.1 03/07/2008 03:33 160 1.5 131 



219+4 
219+4 
210+16 


213+5 
203+6 
200+8 


233+5 
232+4 
213+34 


256+8 
278+9 

237+3 





HD 144432 (calibrator: HIP 018266) 






2117.7 


05/07/2008 02:35 160 1.1 73 


214+3 


220+3 


2123.4 


05/07/2008 02:40 160 1.1 76 


218+4 


230+3 


2236.1 


05/07/2008 03:18 160 1.1 86 


203+25 


229+3 




HD 150193 (calibrator: HIP 090422) 






2117.7 


05/07/2008 05:01 160 1.2 99 


216+5 


192+4 


2123.4 


05/07/2008 05:07 160 1.2 100 


217+3 


197+4 


2236.1 


05/07/2008 05:50 160 1.4 102 


231+27 


206+4 




51 Oph (calibrator: HIP 081266) 






2117.7 


04/07/2008 02:25 80 1.0 -96 


220+3 


221+3 


2123.4 


04/07/2008 02:30 80 1.0 -96 


220+3 


230+3 


2236.1 


04/07/2008 03:06 80 1.0 -101 


215+5 


229+3 




HD 169142 (calibrator: HIP 095347) 






2117.7 


04/07/2008 04:58 160 1.0 63 


208+3 


230+3 


2123.4 


04/07/2008 05:03 160 1.0 66 


219+4 


208+3 


2236.1 


04/07/2008 05:45 160 1.1 81 


216+27 


217+6 




R CrA (calibrator: HIP 095347) 






2117.7 


04/07/2008 08:15 40 1.4 89 


219+5 


230+3 


2123.4 


04/07/2008 08:19 40 1.4 89 


218+3 


208+3 


2236.1 
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217+6 



HD 

03/07/2008 
03/07/2008 
03/07/2008 
HD 190073 

03/07/2008 
03/07/2008 
05/07/2008 



179218 (calibrator: HIP 112029) 
05:29 160 1.3 172 218+3 

05:34 160 1.3 170 222+2 

06:17 160 1.4 157 222+33 

(calibrator: HIP 1 12029 & HIP 090422) 



08: 
08: 
09: 



22 
27 
29 



160 
160 
160 



1.5 
1.5 
1.8 



132 220+6 
131 220+3 
125 223+35 



225+6 
217+3 
220+2 

192+4 
197+3 
206+4 



Notes. *"' Observations from ESO program 079.C-0860(C): (i) observations with the slit aligned north-south; (ii) No suitable PSF standard 
observed at the night of the observations, we used a standard star from our program for the flux calibration. *** Bad weather conditions during the 
observation. 
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Fig.A.l. Observed CRIRES spectrum at the location of the H2 1-0 S(l), 1-0 S(0) and 2-1 S(l) lines in the 13 targets in which IJj 
emission is not detected. The spectra are presented at the rest velocity of the stars. Upper limits to the line fluxes are given in Table 
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